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Abstract. Wireless sensor networks are usually composed of a large number of 
nodes, and with the increasing processing power and power consumption effi¬ 
ciency they are expected to run more complex protocols in the future. These pose 
problems in the field of verification and performance evaluation of wireless net¬ 
works. In this paper, we tailor the mean-field theory as a modeling technique to 
analyze their behavior. We apply this method to the slotted ALOHA protocol, 
and establish results on the long term trends of the protocol within a very large 
network, specially regarding the stability of ALOHA-type protocols. 
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1 Introduction 

The Internet of Things requires the connection of many nodes, therefore we expect that 
in the future we will build Wireless Sensor Networks (WSNs) with a very large num¬ 
ber of nodes. The analysis of such systems dates from 1970’s but still many problems 
have not been solved. In particular, the behavior of such networks under heavy traffic 
loads requires improvements. Typical theoretical studies can only cover a limited num¬ 
ber of aspects, while conclusions from simulations of large networks are often hard to 
generalize. 

In this paper we aim at developing a better understanding of the overall health of the 
network, for instance in terms of the number of nodes that remain in backlog; that is, 
those nodes which have been unsuccessful in delivering their messages and keep trying 
to retransmit the same message. In the past, such analyses have been presented, e.g. [11, 
15], but these studies modeled ALOHA, which is a specific radio protocol that supports 
communication from many nodes to a single central base station. We witness a need to 
model the network protocols in more detail [9]. 

In fact we observed that mean-field theory can be a powerful method to extend the 
approach in [11] and [15]. Two of the most important ways in which mean-field theory 
outperforms the explicit Markov chain analysis is that it does not depend on the size of 
the network, even so, it gives better approximation when the numbers are very large; 
and at the same time with more complex specifications in protocols, the size of the 
semantic model only grows linearly in terms of the number of equations. 

In this paper we show that mean-field theory can reconfirm the results in [8] and 
[15], and aim to show that the mean field approximation method can be suitable to 



model networks. In particular, we focus on an important behavior which occurs in net¬ 
works, called bistability, i.e., the situation where the network may stay for prolonged 
periods of time either in a good state with favorable performance or in an undesirable 
state in which many nodes keep repeating messages, but their transmission are lost be¬ 
cause of mutual interference. That is, the system can converge to two or more different 
steady states, each of which would be a valid solution for the system. To the best of our 
knowledge such behavior has not yet been studied using the mean field approximations 
up until now. 

The mean-field theory originally took shape as a method of approximating complex 
stochastic processes in physics, and was first formally described for communication 
systems in [3]. However, similar heuristics have been proposed before in [10] which 
conformed with steps taken while doing the analysis by means of a continuous time 
Markov chain, in the context of stochastic process algebras. Steps to apply these ap¬ 
proaches to protocols in communication networks have also been taken e.g., in [2]. 
Apart from these more formal and more structured approaches, mean-field analysis has 
also been applied to the analysis of exponential back-off algorithms in communica¬ 
tion protocols [4,5], and to study end-to-end delay and throughput in ad hoc networks 
together with non-equilibrium methods which focus on modeling specific topologies 
[14]. 

For this text, we are going to apply the methods as discussed in [6] to modeling 
wireless sensor networks and do a detailed modeling of phenomena that are observable 
in such networks. The structure of the rest of this paper is as follows: in section 2, 
we give a concise explanation of the essential theory of the mean field approximation 
method, and then in section 3 we first explain interference in networks and then apply it 
to a network running the slotted ALOHA protocol in a simple setting, and consequently 
show its bistability. In section 4, we first describe a way to include broadcasting and 
then consider a slotted ALOHA network with a more complex node behavior. Finally, 
in section 5 we briefly point to the use of vector fields in identifying the conditions in 
which a network shows stability in its behavior. 

2 Mean-field approximation 

In this section, we describe the necessary concepts of the mean field approximation 
method. We consider that every system consists of a number of components. Compo¬ 
nents are formally defined as follows: 

Definition 1. (Component). A component is a quadruple C = (S, Act, — >■, sg), where: 

- S is a nonempty set of states. 

- Act is a set of action names. 

- —>C S X Act X R X S is a transition relation with R = U {_L} called the set 

of transitions rates. 

- Sq G S is the initial state. 

We also write a transition {s,a,r,s') €—)• where a G Act, r G R and s,s' G 
S as s g\ \ more descriptive name for the concept of component would be 



component transition system, since it closely resembles transition systems and other 
familiar automata. 

For a component C, at any point in its computation it can be in a state Y G S which 
is called the state of the component. Now let us consider a system composed of N 
identical components. 

Definition 2. (System). Let Ci = (S, Act, Sq.) where 1 < i < N be identical 
components, except for their initial conditions sq^. A system A4 consisting of these 
components is defined as: 

M = (Ci ||...||CAr). 

Where 11 is a parallel operator. Taking the state of each component Ci to be de¬ 
scribed by a corresponding Yi G Si, we define the state of the system as: 

{Y,,...,Yn). (1) 


Since all components are identical, the mean-field theory allows the modeling of 
system behavior on an abstract level. The overall state of the system is equivalent for 
different permutations of values in the vector (1), and the state of the system can be 
represented in a numerical vector form. 

We take a system model A4 with N components Ci,l < i < N, each with n states. 
The numerical vector form of M., is a vector with n entries which is an alternative 
representation of the state of the system 

=(vr.A-w). 


Here we explicitly mention N in the notation, since in general we take the size of 
a system as a variable. The entry Xj^'^ in vector X^'^ records the number of compo¬ 
nents which are currently in state j. Quick observations are that 0 < Xf^ < N and 


E 


fell, 


.}Af‘=iV. 


Next we describe a normalized vector X^'^ derived from X)A ’ to be the following: 


AN). 




where for each entry X^^'^ = ^ 


(N) 


-. The entries of this vector are also called occupancy 


^3 N 

measures. 

Based on this latter representation of the state of a system we define a class of 
automata called Population Continuous Time Markov chains, which are based on the 
normalized Population Continuous Time Markov Chain models (PCTMCs) with system 
size N, described in [6]. 


Definition 3. (Normalized PCTMC model). Let Ai be a system of N identical com¬ 
ponents C, where the number of states of C is n. A normalized Population Continuous 
Time Markov Chain model for AA is a triple = (X^^^, ), where: 



- is a set of all vectors of the form {x[^'^, ■ ■ ■, Xn^'^), where for all 1 < i < n 

it holds that 0 < x\^^ < 1 and = 1- These vectors are the states of the 

system. 

- X e x(^\ is a vector of dimension n which shows the initial state of the system. 
When the initial state of each component Ci is sq;. d d defined as: 


_ 

— 


yN 


^0 N ^ 

where for b a boolean formula: 

l{b} = 


= Si} 

2^i^l 

5 ■ 

/I : 

. . . , 

ifb true 

|o : 

if b false 


N 


is an indicator function. 

- = {fi,..., fm} is a set of transitions of the form fj = (oj, 

where: 

• Gj € A is the label of the transition. 

• is a vector of dimension n, which shows which portion of each entry in a 
vector X € X^^'^ is going to be used up by this transition. 

• is a vector of dimension n, which shows which values will be added to 
each entry in a vector x € X^^'^ after this transition. 

• : A X X^^^ —)■ M is the rate function of the transition. 

Moreover, we define the state-change vector ~ ) of fj, showing 

the changes in the state of the system due to the transition fj. 

Here we have not demonstrated how is constructed, since it heavily depends 
on the way communications happen in a certain domain. In this paper we will give two 
possible procedures (Transformations 1 and 2) which bridge the gap between compo¬ 
nent transition systems of wireless network nodes and the set of transitions of a 
network of N nodes. 

In the following, the term PCTMC models always refers to normalized PCTMC 
models, unless otherwise stated. Next we give a formal definition of a System of Ordi¬ 
nary Differential Equations (ODEs) for a PCTMC model. 

Definition 4. (System of ODEs). Let x^m ~ , X^^'’ , x\^'' ) be a PCTMC model. 

We consider an n-dimensional vector of functions x(t) = (aii (t),..., Xn{t)) where for 
every t € M>o, x(t) € X^^K The System of Ordinary Differential Equations for the 
PCTMC model x^m defined as: 


^dxift) dXn{t) (tv) (At) 


( 2 ) 


together with the initial condition: 

(ii(0),...,x„(0)) =^0^^ 

where we have f = (a/, r). 


( 3 ) 



The system of ODEs above, together with the initial condition of the PCTMC 
model, can be solved to find a solution (or a number of solutions) which are n curves 
showing the time evolution of the system. 

In this paper we intend to use theorem 4.1 from [6] as a basis for the validity of the 
system of ODEs in approximating the system model. This requires the PCTMC model 
to be density-dependent. Therefore we first present the concept of Lipschitz continuity. 

Definition 5. (Lipschitz Continuity). A function f is called Lipschitz continuous iff there 
is a positive, real constant C, such that for all points x and y over its domain: 

\f{x) - fiy)\ <C-\x-y\. 

Note that any continuous function or any function with a bounded first derivative is 
Lipschitz continuous. 


Definition 6. (Density-dependency). A PCTMC model Xq^'^ ), is 

density-dependent iff the following conditions hold for all x € and f € 

1. for the state-change vector of f, the vector N ■ is independent of N. 

2. there is a Lipschitz continuous and bounded function gj : X^^'> —)■ IR>o such that 
for all X the rate function scales with system size N as: 

rf^\x) = N ■ gf{x). 

Now we present the main theorem (theorem 4.1 from [6]), which states that systems 
of ODEs are approximations for density-dependent PCTMC models when the number 
of components N is sufficiently large. 

Eor every N, let = (S, Q) be continuous-time Markov chains for PCTMC 

models x^m ~ where the set S and the generator matrix Q are 

defined as follows: 

- s = , meaning that and have the same state space. 

- Eor each state € X^^^ and for each transition / = (a, in 

jr(Ar), jf 

x^^'> allows /, meaning that in the vector {x^^'> — there are no 

negative entries, and if (x^^^ -\- € X^^\ then we have the following entry in 

the generator matrix Q: 

Qij = , with i = x^’^^and j = x^^^ -f 

We define the transition probability matrix P^(jv)(f), with entries Pij{f) defining 
the probability of being in state j at time t, if we have been in state i at time 0. Now we 
define the Markov process : f > 0} where = 1, and: 


¥{M, (t) = x^^'^} = pij (t) , where i = Xq^^ and j = x^^^. 



Theorem 1. (Deterministic Approximation for PCTMCs). For values N > lletx^J^ = 
Xq^'^) be density-dependent PCTMCs, andlet (f) : t > 0} denote 

the corresponding Markov processes. Assume that for some point xq € we have 

limjv_>.oo ^ 0 ^^ = io- Moreover, assume that the solution to the system of DDEs for 
X^M ^ (^)- Then for any finite time horizon T < oo, it holds that: 

P{ lim sup I |7W (t) - x(i) 11 = 0} = 1 

Af->oo o<t<T 

The theorem above states that as N grows, the difference between the behavior of 
the explicit model and its approximation x{t) diminishes almost surely. 

3 Modeling interference in Slotted ALOHA networks 

In this section we are first going to discuss interference in wireless networks and how it 
can be modeled and then apply this modeling approach to a Slotted ALOHA network. 

3.1 Interference 

In a WSN an increase in the number of nodes trying to send messages simultaneously 
increases the chances of interference between the signals. In this section, we are going 
to use results from [15] to better express the notion of interference in networks. 

Despite interference, in wireless networks there is always a chance for one of the 
signals to be strong enough to capture a receiver [15]. We consider a set of i senders 
within range of a receiver antenna, and we express the probability of a signal from one 
of the nodes capturing the receiver from a distance r* as: 

q{i) = i ■ V{capture | (z — 1) interfering signals , distance rt} (4) 

Since any of the z messages has a chance of getting through, and these probabilities 
are mutually exclusive, meaning that it is not possible for two or more packets to capture 
the receiver simultaneously (therefore always q{i) < 1). 

We express the scattering pattern of nodes around a receiver as a probability density 
function /(r), henceforth called spatial distribution function, which gives the proba¬ 
bility of a node being at distance r from the receiver. The conditional probability of 
capture for a node becomes [15]: 

roo r ^oo r/ \ 1 (i- 1 ) 

<?(*) = * • / / - 8 — firt)drt (5) 

Jo [Jo l + 2r(r“P . 

Parameter 2 ; describes a power threshold: when facing z — 1 interfering signals, 
the power of a successful message should always exceed the power of the joint z — 1 
interfering signals by z. Parameter /? shows pathloss: the signal gets weaker as distance 
r increases, according to r“^. The term in brackets is the probability of surviving a 
single interfering signal when the power threshold is z and the sender is at distance rt, 



and the enclosing integral averages over different values for rt when we know that it is 
distributed according to spatial distribution function /(r). 

We assume that the spatial distribution function /(r) can have two patterns. A uni¬ 
form f{r) is: 


r 2r for 0 < r < 1, 
0 for r > 1. 


( 6 ) 


Which is the probability of the node being on a circle with radius r around the receiver, 
where r is taken to be a normalized distance. A log-normal fir) is ([15]): 


f{r) 


f} ( log t 


(7) 


Here, cr^ is the spatial logarithmic variance. The log-normal distribution is more 
realistic, in the sense that no node can be infinitely close to the receiver and hence the 
probability of capture when the number of nodes is very high will actually go to zero, 
which is the case in reality. 

Plots of the uniform and log-normal distributions in formula (6) and formula (7) for 
P = A and (7^ = 2 are presented in Fig. 1. 



Log-normal 

Uniform 


Fig. 1: Plot of the uniform and log-normal spatial distributions 


/(r) = 0 and also 

lim^^oo 4 tct = 0. For bounded and positive values of parameter z, the function q{i) 


Lemma 1. Let f{r) be any spatial distribution for which limr_>.oo 
fir? 

in formula (5) is Lipschitz continuous over domain [1, oo). 


Proof Taking the first derivative of formula 5 yields: 
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fir) 
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fir) 


1 -F zrfr~d 
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which we show to be bounded over the domain [1, oo). We establish that the two ad¬ 
dends are bounded, and hence the entire term is bounded. 


- The first addend; the assumption is that 2 ; is positive, therefore for any positive r 
and 

1 + zrf> 1 - ^-5 - < 1 

1 + 2rf r-/5 

Since /(r) is a probability density function, we have f{r)dr = 1. Therefore, 
for any r^: 

r fir) 


lo 1 + zrfr ^ 

r fir) 


dr < 1 


1 (»-i) 


-dr 


^00 

f{rt)drt< / [l]^^-'^'^f(rt)drt = l 

Jo 


lo 

and the first addend is bounded. 

- The second addend; the new term log ( ^r) is unbounded (is - 00 ) 
when; 

r fir) 

Jo 1 + zrf r~^ 

But according to assumptions, since f{r)dr = 1 and limr_).oo fir) = 0 for 
some finite r, f{r) > 0. Therefore —00 happens only when r* —+ 00 . In that case 
it is sufficient to show that; 


0 l+zrfr-^* 
dr = 0 


lim log 


fir) 


-‘-'-'ft \ I a 

rt-j--i-oo Jq l-\-zr^r~l^ 

is zero. We know that; 


dr) ■ f{rt) 


1 


lim 
n-f-i-oo /(rt) 


= +00 


lim 

rt —>+00 


log(/o“ ,JJi-, dr) _ 


^0 l-\-zrf r~^ 

fin) 


00 

+00 


Therefore, there is sufficient ground to use the L’Hopital’s rule: 

poo -Pf{r) 


lim 

—>-+oo 


„ . r^dr f{r,f 


= lim 


/ra r iS*- /'(’■<) 

which is zero since according to the assumptions; 

firt)^ o 

hm = 0 

rt-!-too f{rt) 


and the second addend is also bounded. 


□ 



Corollary 1. Substituting the f {r) functions in formulae (6) and (7) in the definition of 
q{i) in (5) yields Lipschitz continuous functions. 


The interested reader with common skills in calculus could check that the log¬ 
normal spatial distribution in formula (7) satisfies the equation: 


lim 

r—¥(X) 


fir? 

fir) 


= 0 . 


and is hence Lipschitz continuous. As for the uniform spatial distribution in formula 
(6), matters are not as easy, since: 


lim 

r—>-oo 


fjr)^ 

fir) 


- = undefined. 

0 


For the uniform spatial distribution, and for the specific case when /? = 4, we 
propose the following lemma. 


Lemma 2. The function: 




2r 


n (*-l) 


+ zrfr~'^ 


dr 


2rtdrt 


= i ■ 


to L 


1 — s/zr'l arctan( „ ) 


2rtdrt 


is Lipschitz continuous. 

Proof. The first derivative of q{i) is: 


q'{i) = 


1 


-| (*-i) 


1 — f zr^ arctan( „ ) 


2rtdrt 


+ i ■ 


1 — s/zrf arctanf ^ ), 

W^rfl 


(*-i) 


log (1 — fzr^ arctan{—=-^))2rtdrt 

y^r^ 


For the first line we know 0 < arctan(a;) < f , a: G K>o; therefore for i > 1 the 
first factor is bounded, and by integrating over the bounded domain [0,1], it will still be 
bounded in value. 

For the second line, we again have the same term multiplied by: 


1 


log (l — -s/ir? arctanl ^ 
yfzrt 


which is independent of i. This term is unbounded (evaluates to —oo) only when: 


1 


(l — \/zr1 arctan( ^ 2 )) “ 



1 


arctanf „ ) = 

Vzrt 


Vzrf 



= 0 . 


But this is not the case here since z is bounded and 0 < < 1. After integrating over 

the bounded domain [0,1] and multiplication by i the term is still bounded. □ 


In Fig. 2, values numerically derived for formula (5) are shown when the number 
of signals is between 1 to 10, for /3 = 4, z = 10 and = 2. For uniformly distributed 
nodes (the points are shown with ‘+’), the curve will stabilize over a fixed value (~ 
0.2) for large numbers. But when using log-normal spatial distribution to calculate the 
capture probability q{i) (the points are shown with ‘o’), the curve goes to zero. 



Fig. 2: The probability q{i) for 0 < i < 10. The two set of points show two different patterns: q{i) 
with log-normal spatial distribution (o) tends to 0 for very large i, while with uniform distribution 
(-I-) it converges to a fixed value (~ 0.2). The lines show a possible interpolating function as a 
continuous approximation. 


In the rest of this paper, function q{i) is used over the domain [0, oo). So we assume 
that the value of q{i) in the range [0,1) is defined as: q{i) = i, since we know that for 
a single transmitting node (or less) there is no interference in the network. This also 
assures continuity, since ( 7 ( 1 ) = 1 in both functions. 

The following summarizes the results of this section in terms of the definition of the 
set T in PCTMC models for WSNs. 


Transformation 1. For component Node = (S', Act, —>■, Sq), assume that the transi¬ 
tions in the set —> are restricted to the following forms: 


Si 


{capture send] 


> Sj be “successful send” transitions. 


(failure,r send) 


> Sfc be “failed send” transitions, due to interference. 


{cx,ra) 


- Si 


■> Sm €—)■ be any other transition. 



Let Network = {Node[^'^ 11 • • • 11 Node \!^'^) be a system made of N identical compo¬ 
nents of type Node. The set of the associated PCTMC model xj^tlork ~ 

defined for x, where x € consists of the following tuples: 

ffafiL^ = (failure, Pseud -{N •x^- q{N • Xi))) 

f capture.^ = {capture, Psend ■ q{N ■ xf) 

/<^f^ =ia, r^-N-xi) 


where is a unit vector with value 1 at position i. 

Lemma 3. Let N > 1 and function q{i) be a Lipschitz continuous function over do¬ 
main i G (0,oo) and right-continuous at point i = 0. The PCTMC model xl^tlork ~ 
{X^^'l, is density-dependent. 

Proof. In order to prove density-dependency for X^^etwork check the conditions in 
Definition 6 : 

1. The vectors: 


= (cj — ef ), 

,Ck C?' 

N-U(n) = N ■ {—-—) 

•'failure^ iV iV 


— (Cfc Ci) , 

— (^m 




are independent of N. 

2. As for the continuity criteria, we consider the gf functions below to be functions of 
the entries of the vector X^^'i and prove their Lipschitz continuity: 

dgpN) 

- The function g ^{n) = Va ■ xi has the first derivative = r^, where Tq 

is constant and therefore continuous and bounded. This means that g an) is 

Jai 

Lipschitz continuous. 

- For the function g an) = Vsend ■ jj ■ q{X ■ x,), the term rgend ■ jj would be a 

d failure^ 

constant for a predetermined value of iV > 1. As for the second part, we know 
that Q < Xi < 1. But then 0 < N ■ xt < N since > 1 according to the 
assumptions. Again, based on the assumptions we know that q{i) is Lipschitz 
continuous over domain i G [0, oo). Therefore the term Pan) is Lipschitz 

J failure j 

continuous. 



- The function (/,(») = Tsend ' {xi — ^ ■ q{N ■ Xi)), has two factors, the first 

J capture j 

of which is 

Xsend ’ Xi 

and is Lipschitz continuous like the first case. And the second term 


Xsend 

N 


■ q{N ■ Xi) 


is Lipschitz continuous following the second case. 


□ 


3.2 Slotted ALOHA with a single receiver 

In this part, we will consider a PCTMC model, built according to Transformation 1, 
where every Node component in the system runs the Slotted ALOHA protocol. We 
consider the same scenario as in [15] which consists of a number of senders scattered 
around a single antenna. The aim is to use the mean field approximation method to 
observe the bistable behavior of this specific ALOHA network. 

The Slotted ALOHA protocol we consider here is expressed by the following set of 
rules [1,12]: 

- Whenever there is data to send, send it at the start of the next time-slot. 

- If the message could not be delivered due to interference, retry sending the message. 

To this we also add the following restriction: 

- While sending and retrying, do not generate new messages. 

A node’s behavior is presented in Fig. 3. A node does its internal processing in state 
(O), and generates a new message with rate Vo- Next, the node enters a state where 
it transmits the message (T). While sending with rate rggnd. if other nodes are also 
transmitting messages simultaneously, the signals will interfere. In case the message 
cannot be delivered, a node enters the backlog state (R), where it tries to retransmit the 
message after some time (r^). 


(generate, To) (failure, Vsend) 



(capture, Tsend) (reseud, r^) 


Fig. 3: An ALOHA node’s component transition system 


Now we proceed to generate the set of differential equations associated with the 
slotted ALOHA network. Following the theory presented in section 2, the number of 
processes in state s G {0,T,R} at time t is expressed by the function Xs{t), and 
following Transformation 1 we derive the set of equations in Table 1. 



Table 1: System of ODEs for the Slotted ALOHA model 


■ Xo{t) + ^ • fsend • q{N ■ XT{t)) 

= -»'send • (XT(t) - ^ • q{N ■ XT{t))) - ^ • r^end ' q{N ■ XT{t)) 
+ To ■ Xo{t) + Tr ■ XR{t) 

= -rr ■ XR{t) + Tsend ' (*T (f) “ ^ ' q{N ■ XT(t))) 


Following [15] we take = 0.08, Xq = 0.0055 and rgend = 1; since a message 
transmission always takes one timeslot. Specifying the initial condition as (i;o(0) = 
a:T(0) = 0, Xi?(0) = 0), and solving the equations numerically (because the complexity 
of function q{i) does not allow explicit solutions) we have the curves in Fig. 4a. Here 
the curves show the number of nodes in each state, which change over time. After some 
changes, they stabilize over an equilibrium point, which is the fixpoint of the system of 
ODEs. The plot shows a network with a good behavior, in which once in a while a new 
message is generated and is almost always successfully delivered to the receiver in the 
first try {xR{t) ^ 0). 

However, under different initial conditions, namely (i(3(0 ) = 0, xt( 0) = 1, a;ij(0) = 
0), the solution is the curves in Fig. 4b. This plot shows that if the system starts in a state 
where everybody is trying to transmit a message, then it will be trapped in a state with 
a very low throughput where a constant number of nodes are always trying to deliver 
their messages and saturate the media. 

In Fig. 4c, we see that in a less realistic case of traffic in a uniform spatially dis¬ 
tributed network, the nodes tend to operate efficiently after some time despite the initial 
conditions. 

The bistability that we see in case of the log-normal spatial distribution is a common 
property of many ALOHA-type networks. The results match those presented in [15], 
and other observations on how ALOHA networks behave in real world [13]. 

4 Modeling local broadcast in Slotted ALOHA networks 

In this section we are first going to discuss local broadcast in wireless networks and 
how it can be modeled, and then apply the modeling approach to a network running a 
simple neighborhood discovery protocol. 

4.1 Local broadcast 

In essence, wireless network communication traffic consists of a number of broadcasts, 
in which every node in the vicinity of a sender is capable of hearing the message; so a 
message transmission typically involves one sender and several receivers. 

The number of receiver parties depends on the technology of the radio modules, the 
properties of the media, the network topology, etc. Here, we consider a single measure p 
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Fig. 4: Numerical solutions for the slotted ALOHA network 


to represent all these properties. The probability p is the fraction of nodes that are close 
enough to the transmitter to be able to successfully receive the message, and usually 
has the form p = where d is some constant number of nodes which are located in a 
certain neighborhood. 

For this communication, two terms appear in ODEs; one for the receivers and one 
for senders. For the receiver, we start with the term used in [7] to describe communica¬ 
tion when there is only one receiver. We call this M^, which has the general form; 

= ra • Ai, • l{Nr > 0}. (8) 

Where is the rate of send action a, Ns stands for the total number of processes in 
states in which a send action a is possible, Nr is the total number of processes in states 
in which a receive action is possible, and the indicator function l{Afr > 0} shows the 
fact that communication is not possible when there are no receivers. 

Instead, in broadcasting we reintroduce the term Ma in (8) as term 11 a. 

Ua = ra ■ Ns ■ {p ■ Nr). (9) 

since at each point in time we know that a portion of all receivers that are within range 
(expressed by parameter p) are capable of receiving the message. 

TTq, in equation (9) can be derived by rewriting the term Mq, = ra-Ns-l{Nr > 0} in 
a setting with multiple receivers. In order to show this, we first describe ¥{M = i\Nr} 













as the probability of i successful receives, when there are a total of Nj. receivers: 


P{M = i\Nr} = (10) 

Because we know that any subset of the possible receivers are capable of participating 
in the wireless communication, we extend term (8) and derive: 

N 

77„ = ^ P{M = -i-ra-Ns- l{Nr > i}. 

Where N is the total number of nodes in the network. Replacing the term P{M = 
z|AV} with the right hand side of (10) we have; 

Nr / \ 

na= ^ jp*(l• i • ra • iVs 
i=o \ ^ 

Nr 

= r^■N,■{J2^rnp\l-pr^-^). 

i =0 ^ ' 


The series is the expected value of a binomial distribution; therefore the receive 
term is simply: 

Ua = ra ■ Ns ■ {p ■ Nr). 

Next, we see how results presented regarding interfering signals can be combined 
with local broadcasting and help us derive the term Ra- Term (9) can be interpreted 
differently, as a portion of senders (p ■ Nr) which are within range of each receiver, and 
try to capture it: 

Ta- (p- Ns) ■ Nr. 

And then we apply the function q{i) to the total number of received messages at 
each receiver’s site: 

Ra = Ta ■ qip ■ Ns) ■ Nr. (11) 

This is the term we intend to use for receivers. As for the senders, since a node which 
is sending a message type a does not depend on the status of receivers, we have the 
following simple format: 

Ta = ra-Ns. (12) 

This means that when transmitting, the sender does not depend on the number of re¬ 
ceivers in its vicinity. 

In the following, all send or receive actions have the form send{m) or receive{m) 
and are accompanied by a type m G M, where M is a set of message types. We define 
the set Ira to contain send actions that are interfering with send actions with message 
type TO, where at least sendim) G Im- 

Transformation 2. For component Node = (S', Act, —>■, sq), assume that the transi¬ 
tions in the set —> are restricted to the following forms: 



(send(m),r„nd) ^ ^ 

Si - > Sj be send transitions of message type m, 

(receive(m),±) ... . . „ 

Sfc- > Si €—>■ be associated receive transitions of message type m. 


- Sm —^- y So be any other transition. 

We forbid the component Node to be able to both send and receive a message type m 
when in a state Si G S. 

Let Network = {Node^.^^ II • • • II Node^j^'^) be a system made of N identical Node 
components. Let be the associated PCTMC model. 

Consider C (^) to be the portion of nodes that are in a state in which they are 
capable of doing the actions of set A, where for G it is defined as follows: 


C{A) = 


E 


AN) 


iG{a|3sb.(sa" 


>SbAaGA)} 


The set consists of the following tuples (again for x defined above): 


AN) 

J send(m)i 


e' e' 

= {send(m), N-rsend-xf) 


fLZie(m)^ = {receive(m), 


N- 


^ send 


C{{send{m)}) 

Cilm) 


q{N ■ p ■ C{{send{m)})) ■ Xk) 


, Cm Co 

= aT’ 


fiV = -Xr,- Xr. 

where Ci is a unit vector with value 1 at position i. 


It is worth noting that the term N ■ p = N ■ = d, and therefore f 


(N) 

f receive(m)^ 


does 


not depend on N. 

Lemma 4. Let N > 1 and function q{i) to be a Lipschitz continuous function over 
domain i G (0, oo) and right-continuous at point i = 0. The PCTMC model xAlvork ~ 
(^X^N) ^ jr(tv) ^ Xq^'^) is density-dependent. 

Proof. We check the conditions in Definition 6; 

1. The vectors: 


N ■ v .(N) 


•> send(Tn.).j 


N • uj{N) 


= 7V.(^--) 

— Ci , 

t Ar ]xf> 


receive(m)^ N 

= Cl — Ck , 


N ■ vjiN) 


= N-A-^AA) 
N ’ 


- 6 o Cr] 


are independent of N. 



2. As for the continuity criteria, we consider the gf functions below to be functions 
over vectors x € and prove their Lipschitz continuity: 

- (N) ^3 JN) 

- The function = Tq • Xm has the first derivative —= Tq, where Tq 
is continuous and bounded. Therefore p^cjv) is Lipschitz continuous. 

- For the function 

C({send(m)}) . ,, .. 

gAN) = Tsend ' - , - q[N ■ p ■ C {{send{m) })) ■ Xi 

Jreceive^ C(Im) 

We first make an observation. We know that 

ieB 

where B C n}. Since sendim) G Im we always have: 

C{{send{m))) < C(Im) 

Next, we show that for two n-dimensional vectors "af and l/ in X^^\ there is 
a constant real number c for which: 

\gAN) {'^)-gf(N) (V)! < c-- Vl 

*'rei 


V, 


holds, where |lzf — ^| should be interpreted as the distance between the end¬ 
points of the two vectors. Returning to the observation that we made at the start 
of the proof, and the fact that q(i) is a probability, we have: 


C{{send{m)}) 


< 1 


and. 

Therefore: 


C{Im) 

q{N ■ p ■ C{{send{m)})) < 1 

C({send{m)}) ,, 

• q{N ■ p ■ C{{send{m)})) < 1 


C{Im) 

and multiplying the two sides by Vgend ■ Xi we have: 
C{{send{m)}) 


^send ’ 


C{Irr 


q{N ■ p ■ C{{send{m)})) < Xaend • Xi 


So in order to establish the main result it suffices to show that: 

l^send ’ i send ‘ i | ^ C * \~^ \ 

for any positive Vgend we have: 

I'^send ‘ i send ‘ i\ — ‘f'send * \~^i i\ 



and since by definition the distance must have the following property: 

we have: 

t* send ‘ I i "l/ i\ — send * I I 

by transitivity of the inequality relation < we have shown that at least for c = 

^ send' 

IgAN) (af)-p,(N) (V)l ^ ^send l^-Vl 

d receive^ *' receive^ 

The inequality holds and the function g an) is Lipschitz continuous. 

•'receive^ 

- The function g an) = r^end • Xi, is Lipschitz continuous following the proof 

•' send^ 

for the first case. 
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4.2 Neighborhood discovery protocol 

In this part we will study a slightly different version of the neighborhood discovery 
protocol in an ALOHA network. The discovery works as follows: every once in a while 
each node broadcasts a HELLO message to advertise its presence in the network. All 
the neighbors hearing this will respond with an acknowledgement. The sender follows 
a passive acknowledgement model, and upon receiving an acknowledgement which 
ensures its discovery by at least one neighbor, proceeds with its internal processing. 
Such protocols are essential building blocks of many algorithms such as routing in 
wireless ad hoc networks. 

We consider every node to run an identical implementation of the neighborhood 
discovery protocol. An abstract transition system for every node’s behavior is given in 
Fig. 5. 



Fig. 5: Component transition system of ALOHA nodes in the neighborhood discovery protocol. 


In this transition system multiple assumptions and considerations are made. Re¬ 
ceiving a message is allowed at any point in time when a node is not busy sending a 



message. Therefore when a node is doing its internal process, or when it is waiting for 
an acknowledgement, it responds to any message that is received in the mean-time. 

Timeout, which is an internal action, is taken to have a fixed duration {T) and the 
rate of the timeout is defined accordingly as (rtimeout = ^)- Nodes do other processing in 
between sending messages, which also happens after a fixed —-— time. A node which 
has a good performance spends most of its time processing in state 0, instead of trying 
to send a message. 

Based on the transition system in Fig. 5, and also that Imsg = ack} and 

lack = {msg, ack}, we derive the system of ODEs according to the recipe provided in 
Transformation 2. The result is presented in Table 2. 


Table 2: System of ODEs derived from the transition system in Fig. 5. Vack and rmsg are only 
used for clarity and are otherwise equal to Vsend- 



We take the parameters (vsend = 100, rprocess = 1, rtimeout = 30) and p = 0.05 
for the connectivity of links in the network. We take N = 500, and therefore we have 
the number of nodes in a neighborhood d = • p = 25. By solving the system of 




equations for the initial conditions (io(0) = 1) and then for (i 3 ( 0 ) = 1), we see a 
bistable behavior, as can be seen in Fig. 6a and 6b. 



(a) Solution with = 1 (b) Solution with is = 1 

Fig. 6: The bistable behavior in an ALOHA network with multiple receivers 


5 Stability and Vector Field Analysis 


A common method to study stability in Markov counting processes has been the calcu¬ 
lation of drifts. However, as the model grows, observing drifts in differential equations 
is not straightforward. Therefore we use vector fields to study the equilibrium points. 

For equations in Table 1, a corresponding two dimensional vector-field has been 
given in Fig 7a. Here, since we know that xq + Xt + in = 1, only the triangular area 
at the bottom of the plane is filled, where io + xr < 1. The value of xt at each point 
in this area is implicitly defined as xt = ^ — xq — xr. The vectors in this figure are of 
two colors, and show the tendency of the system starting from that point (as the initial 
condition) to go in any two directions; the red part consists of points which go to the 
first fixpoint, and the blue part consists of points which go to the second one. 

For equations in Table 2, Fig. 7b shows a similar pattern. Here, since there are 6 
varying parameters Xq through x^, we only present two values for each point, one for 
parameter io which is the original state and the other for x^, the state in which a node 
retries sending a messages. Also we take io + ^3 < therefore the rest of the 

parameters may take positive values which satisfies the equation ii -F i2 -I- i4 -F io = 
1 — io + is. For each point. Fig. 7b shows the average tendency to go to either of the 
two fixpoints, where white areas go to the first fixpoint and black areas go to the second 
one. For areas in which it is possible to go to both fixpoints for various values of the 4 
absent parameters, the color will turn out to be different shades of gray, depending on 
the proportion of cases which go to each fixpoint. 

In most systems, avoiding a bistable behavior improves the predictability of the 
system. However, when bistability is inevitable, a system’s conditions can be monitored 
for signs that make it prone to slipping into a state with low performance, e.g. when the 
system in Fig. 7b is in the gray area. 













Xo ^ 

(a) Tendency of each {xo,Xr) pairs to go to (b) Tendency of each (xq, £3) pairs to go to 
either of the fixpoints (disks). either of the fixpoints. 

Fig. 7: Visualization of bistable behavior in ALOHA 


6 Conclusion 

In this paper, we took crucial steps in modeling large wireless networks. We used in¬ 
sights from the field of wireless communication to describe communication models, and 
customize the mean-field theory. In this way we are able to reason about networks which 
are immensely large, and avoid common limitations in analyzing traditional models. 

The reasons for the correctness of this modeling approach is twofold. First, we 
proposed semantics (Transformation 1 and 2) which fully conformed with theories 
that guarantee effective approximation of Markov chain models with systems of ODEs 
(Lemma 3 and 4). Second, through examining ALOHA networks we were able to wit¬ 
ness the same bistability phenomena that are commonly associated with them in prac¬ 
tice. 

Linally, we demonstrated the modeling capabilities of this approach with a sim¬ 
ple neighborhood discovery protocol. Moreover, we focused on systems with multiple 
equilibrium points and used vector-fields as a way to visualize the dynamic behavior of 
the system. 
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